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(N ; We discuss the possibility of using electroproduction of J/ip as a probe of gluon Sivers function 

by measuring single spin asymmetry (SSA) in experiments with transversely polarized protons and 
electron beams. We estimate SSA for JLab, HERMES, COMPASS and eRHIC energies using the 
color evaporation model of charmonium production and find asymmetry up to 25 % for certain 
choices of model parameters which have been used earlier for estimating SSA in the SIDIS and 
Drell-Yan processes. 
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I. INTRODUCTION 



In recent years, there has been a lot of interest in investigations of transverse single spin asymmetries (SSA's) in 
high energy QCD processes as they provide information about spin structure of proton. These asymmetries arise in 
scattering of a transversely polarized proton off an unpolarized proton if the scattering cross section depends on the 
\ direction of polarization. 
\Q . The single spin asymmetry (SSA) for the inclusive process + B — > C + X is defined as 

O 1 



A N = -7-^-J-T (1) 



da^ + dcr-l- 

£S| . where dcA^ denotes the cross section for scattering of a transversely polarized hadron A off an unpolarized hadron 
B, with A upwards (downwards) transversely polarized with respect to the production plane. Large SSA's have been 
measured in pion production at Fermilab [l| as well as at BNL-RHIC in pp^ collisions [2|. SSA's have also been 
observed by the HERMES [1] and COMPASS 0] collaborations, in polarized semi-inclusive deep inelastic scattering. 
The magnitude of the observed asymmetries have been found to be larger than what is predicted by pQCD Q . 

It was first proposed by Sivers Q that it may be possible to explain this asymmetry by allowing a correlation 
between the transverse momentum of the quark and the polarization of parent hadron. This approach requires a 
pQCD factorization scheme that includes the spin and intrinsic transverse momentum effects. With the inclusion of 
fc_i_ dependence in parton distribution functions (pdf's) and fragmentation functions (ff's) Q, one is led to a generalized 
factorization formula called transverse momentum dependent TMD factorization [f| H| . TMD factorization in some 
processes has been proved at leading twist and leading order [§] and has been argued to hold at all orders. 

The inclusion of the effect of transverse momentum of partons in pdf's and ff's leads to a new class of parton 
distributions which are obtained as extensions of usual collinear pdf's and include the effects of spin and transverse 
momentum of the partons. One of these functions is the Sivers function, which describes the probability of finding 
an unpolarized parton inside a transversely polarized hadron. The coupling of the transverse momentum of the 
unpolarized quarks and gluons to the nucleon spin is in fact related to their orbital angular momentum. Thus the 
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Sivers asymmetry gives access to the orbital angular momentum of the partons. The number density of partons inside 
a proton with transverse polarization S and momentum P is parameterized as [To| 

fa/p^(Xa,k± a ,S) = fa/p(Xa,k± a ) + - A N f a/pt (x a , fc_L a ) S ■ (P X fe± ) (2) 

where k± a is the transverse momentum of the parton, x a is the longitudinal momentum fraction of parton, 
fa/p{ x jk± a ) is the unpolarized parton distribution and A N f a / p t (x, k± a ) is the Sivers function. In this work, we 
propose charmonium production as a probe to investigate the Sivers function and as a first step, estimate SSA in 
photoproduction of charmonium in scattering of electrons off transversely polarized protons. 

One of the difficulties in getting information about the spin and transverse momentum dependent pdf's and frag- 
mentation functions is that very often two or more of these functions contribute to the same physical observable 
making it difficult to estimate each single one separately. It has been shown how properly defined SSA's in the 
Drell-Yan process depend only on the quark Sivers function and the unpolarized quark distributions [lol ] . The studies 
of Anselmino etal show that the magnitude of the Sivers asymmetry in the Drell-Yan (DY) process for forthcoming 



experiments at RHIC, COMPASS, J-PARC, PAX, PANDA, NICA and SPASCHARM is large A study of the 
Sivers effect for pion and kaon production in semi-inclusive deep inelastic scattering (SIDIS) processes has been per- 
formed and estimates have been given for experiments at COMPASS and JLab [12l4l4| . It has been proposed to probe 
the gluon Sivers function by looking at back-to-back correlations in p^p interactions at RHIC . Another process 
that has been suggested as a probe to access gluon Sivers function is p^p — > DX as the SSA in this case arises due to 
the gluon Sivers function alone mainly in the intermediate rapidity region (l6| . 

In this paper, we investigate feasibility of using charmonium production to obtain information about the Sivers 
function. Charmonium production has been known to be a sensitive tool to study QCD for bound states of heavy 
quark-antiquark systems. Recently, there has been some discussion about the possibility to use fixed target experi- 
ments at LHC for charmonium production with the aim of investigating the quarkonium production mechanism [l7| . 
Photoproduction and electroproduction of charmonium near threshold are expected to throw some light on this mech- 
anism as well as on hadron structure [l8j . In fact the connection between charm, charmonium production and the 
gluon densities has been explored since long for protons-polarized and unpolarized-, nuclei and photons [l9j . Here, 
we study asymmetry in photoproduction (i.e. low virtuality electroproduction) of charmonium in scattering off po- 
larized protons. At LO, this receives contribution only from a single partonic subprocess 73 — > cc . Hence, SSA in 
e + p^ — > e + J/ip + X , if observed, can be used as a clean probe of gluon Sivers function. In addition, charmonium 
production mechanism can also have implications for this SSA and therefore, its study can help probe the production 
mechanism for charmonium. 

There are three models for charmonium production. In the color singlet model 20] the cross section for charmonium 
production is factorized into a short distance part for cc pair production calculable in perturbation theory and a 
nonperturbative matrix element for the formation of a bound state, which is produced in a color singlet state. In 
color evaporation model, first proposed by Halzen and Matsuda [2l| and Fritsch [22j a statistical treatment of color 
is made and the probability of finding a specific quarkonium state is assumed to be independent of the color of heavy 
quark pair. In later versions of this model it has been found that the data are better fitted if a phenomenological 
factor is included in the differential cross section formula, which depends on a Gaussian distribution of the transverse 
momentum of the charmonium 23]. A more recent model of charmonium production is the color octet model pij ]. 
This is based on a factorization approach in nonrelativistic QCD, and it allows cc pairs to be produced in color octet 
states. Here again, one requires knowledge of the nonperturbative color octet matrix elements, which are determined 
through fits to the data on charmonium production. Inclusion of nonzero intrinsic transverse momentum for the 
colliding partons can help us to understand the discrepancy between these matrix elements determined from the 
hadroproduction (Tevatron) data and the leptoproduction (HERA) data [2|| l26j . Since the nonrelativistic QCD 
calculations are done in collinear approximation, it is not surprising that the effects of nonzero transverse momentum 
of the colliding partons can be large. In this work, we have chosen to work with the color evaporation model (CEM) 
as its simplicity makes it suitable for an initial study of SSA in the charmonium production. 

It has been proposed that the SSA's at leading twist in pQCD, arise from the final state interactions between the 
outgoing quark and the target spectator system [271 ]. In case of the Sivers asymmetry, the initial and/or final state 
interactions can generate a nonzero phase in the amplitude, which through the naive time reversal odd (T-odd) Sivers 
function then gives rise to the SSA. This interaction is part of the gauge link present in the TMD functions and 
depends on the specific process under consideration. Thus it introduces a process dependence, in particular, in the 
T-odd distribution functions [28[ . In fact the direction of the link is opposite in the Drell-Yan process as compared to 
SIDIS. For a generic hadronic process, it can be completely different from the above two processes. The generation of 
SSA in J /ifi production in ep^ and pp^ was studied in [29| . It was found therein that existence of nonzero SSA depends 
on the production mechanism of J/ip, namely in ep^ processes nonzero SSA is expected only if the charmonium is 
produced in the color octet state and for pp^ it should be produced in the color singlet state. However no numerical 
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estimate of the asymmetry was given. 

We provide a first estimate of the SSA in the low Q 2 electroproduction (photoproduction) of J ftp at leading order, 
using a factorized formula along with models of Sivers function used in the literature. As a first step in our investigation 
of SSA's in charmonium production, we have used the color evaporation model for charmonium production. According 
to CEM, the cross section for charmonium production is proportional to the rate of production of cc pair integrated 
over the mass range 2m c to Ivnn [30l | 

1 f 2mD ... da c . 



a = - \ dM ce — (3) 

9j 2ma cc dM cB ( > 

where m c is the charm quark mass and 2mo is the DD threshold. 

At leading order (LO) there are no ff's involved and the only contribution to asymmetry comes from the Sivers 
function. Therefore, we can use this observable to extract information about it. 

Here, we have used Weizsacker- Williams equivalent photon approximation to calculate cross section for the process 
e + p^ — > e + J/ip + X . The underlying partonic process at LO is jg — » cc and therefore, the only k± dependent 
pdf appearing is the gluon Sivers function. For a complete calculation of photoproduction of J ftp one has to consider 
higher order contributions and also the resolved photon contributions [23j . 

To assess numerical estimates, we have used and compared two models of the quark Sivers function obtained from 
SIDIS data and have estimated the magnitude of asymmetry at JLab, HERMES, COMPASS and eRHIC center of 
mass energies. We predict nonzero asymmetry in both rapidity and qt distribution for the parameter set fitted from 
experimental data [3l|. 



II. FORMALISM FOR ASYMMETRY IN J/ip PRODUCTION 

A. Color Evaporation Model 

We consider the LO parton model cross section for low virtuality electroproduction (photoproduction) of J ftp within 
color evaporation model. According to CEM, the cross section for charmonium production is given by Eq. ([3]), where 

d(7 C c 

— is calculable perturbatively, M cS being the invariant mass of the cc pair. The differential cross section for the 

dM cS 

7 + p — > J/ip + X is given by 

where f g / p (x) is the gluon distribution in the proton. 

Using the Weizsacker- Williams approximation [32|, [33[ , one can convolute the cross section given by Eq. (j4| with a 
photon flux factor to obtain the electroproduction cross section for e + p — > e + J/ip + X 

= Jdy f lh {y) (5) 

where y is the energy fraction of electron carried by the photon and / 7 / e (y) is the distribution function of the photon 
in the electron given by 34| , 
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f ( ^ < V + (l-j/) 2 f, E V 

f 1 /e{y,E) = -{ In 

tt y \ m 2 , 
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Thus, the cross section for electroproduction of J/tp using WW approximation is given by, 



4™ n , d&19 - 



4 m- 



a e P ^ + J N , + X = I ^ I dydx f i/M fg/p{x) (7) 
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B. Single spin asymmetry in J/tjj production 



To calculate SSA in scattering of electrons off a polarized proton target, we assume a generalization of CEM 
expression by taking into account the transverse momentum dependence of the Weizsacker- Williams function and 
gluon distribution function 

-jj^ = / dx 7 dx g [d 2 k^d 2 kt g ] fg/pr(xg,k ±g )f 7/e (x 1: k^)—jj^- (8) 

where M 2 = M 2 S . We have not written the scale dependence of the quantities on the right hand side explicitly. As 
mentioned earlier, a generalization of factorization formula involving TMD pdf 's and ff 's leads to nonzero SSA in DY, 
SIDIS and other processes. Therefore, we expect that inclusion of transverse momentum dependence in WW function 
and generalization of the CEM expression might be a valid approach to estimate SSA in charmonium production. 

SSA is defined as in Eq. (1), where in our case da™ are single transverse spin-dependent cross sections for 
e + — > e + J/ip + X and e + p^ — > e + J/tp + X, respectively. 

The difference in do^ and da^ is parameterized in terms of the gluon Sivers function 



da 1 ' — do^ = J dx~f dx g d 2 k±~ t d 2 k± g A N f g / p t (x g , k± g ) / 7 / e (x 7 , k± 7 ) do 19 ^ cc (9) 
where do is the elementary cross section for the process 73 — > cc given by 

1 d 3 p c d 3 p E 1 



oV 9 ^ c5 (M 2 ) = ± -J± S((q - Vc f - ml) I M 7ff ^ cc - 1 2 (12) 



d& = Ys 2E c 2E; J2^ S (P7+P 9 -Pc-P 5 )|M 7 ^ CC -| . (10) 
We rewrite t^Jt" — d 4 pc S(p 5 2 — m 2 ) and change the variable to q = p c + Pc so that 

f^ = d 4 qS((q~ Pc ) 2 -m 2 c ). (11) 

Now using the expression for total partonic cross section 

1 f d 3 p c 1 

and changing the variables from qo and q^ to M 2 and rapidity y so that 

dM 2 dy = 2dq dq L (13) 

We finally obtain 

dydM 2 d 2 q T ~ dydM 2 d 2 q T 

xf l/e {x 11 k^)5\p g +p 1 - q) &2 9 ^ C5 (M 2 ). (14) 

where the partonic cross section is given by [35| 



= \ J[dx 7 d 2 kj_ 7 dx g d 2 kj_ g ]A N f g/p r{x g ,kj_ g ) 



Am 2 

Here, 7 = — ^ and M 2 = s. 

A N f g / p t (x g , k± g ) is related to the gluon Sivers function A N f g i v t (x, k± g ) by 

A N f g/pt (Xg,k ± g) = A N f g/p ,(Xg,k ± )S-(P X k ± g ) . (16) 

Similarly the total cross section is given by 

dV dM 



j ' [dXjd"k±ydXgd kj_g]fg/p(Xg,kj_g) 



dydM 2 d 2 q T dydM 2 d 2 q T 

xf l/e (x 7 ,k ±7 )5 4 (p g + Pl - q)o^ c5 (M 2 ). (17) 
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It is worth mentioning at this point that following the procedure used by Anselmino etal in case of Drell-Yan process, 
we have been able to write the distribution in M 2 , qx and y of the produced J/ip in terms of unpolarized total partonic 
cross section. 

Rewriting the four momentum conservation S function as 

^{Pg +Pi - l) = S{Eg+E 7 -q )5(p Zg +p z _ i -q L )5 2 (k ±1 + k ±g -q T ) 

2 J ] ' y ) s( Xj -^e-As 2 (k Xy + k ±g -q T ) (18) 



one can perform the x 1 and x g integrations to obtain 



and 



dydM 2 d 2 q T dydM 2 d 2 q T s 



d 4 aT d 4 ^ 



- / [d 2 k ±1 d 2 k ±g }A N f g/pt (x g 

x 5 2 (k ±1 + k ±g - q T )a^ c£ (M 2 ) (19) 



dydM 2 d 2 q T dydM 2 d 2 q T s 



with 



J \d k±^d k± g ]f g /p(x g ,k± g )fy/ e (x.y,k±^) 

x 5 2 {k^ + k ±g - q T )a^ c5 (M 2 ) (20) 

M ± 

x g,y = ~F e y - (21) 

d 3 a^ 

Integrating Eqs. (fT9l) and (1201) over M 2 as prescribed by CEM, we obtain the difference and sum of - — = — and 

dyd z q T 

d 3 a l 

- for J/ip production. 

dyd 2 q T ' 

We follow the convention in Ref. [36| and dehne the Sivers asymmetry integrated over the azimuthal angle of J/tp 
with a weight factor sin(</> gT — (f>s) ■ 



i(<V-</>s) _ J d(j> q T [da f - da±]sm(<j> qT - <j) S ) 
J d<fi qT [da^ + da^] 



where do^ is differential cross section in qx or y variable and 4> qT and 05 are the azimuthal angles of the J/ip and 
proton spin respectively. To evaluate asymmetry in y distribution, we will substitute 



da^ - da l = f dcj) qT f q T dq T f ° [dM 2 ] f [d 2 k± g ]A N f g/pt {x g ,k± g ) 

J J J Ami * 

x f 1 /e(x 1 ,q T - k ±g ) o- (M 2 ) sin(<p qT - <p s ) (23) 

and 

da^ + da 1 = 2 J d<j) qT J q T dq T J ^ [dM 2 ] J [d 2 k ±g ]f g/p {x g ,k ±g ) 

x f 7 /e{x-y,q T - k ±g )a Q {M 2 ). (24) 



Thus at LO, the SSA depends on Weizsacker- Williams function, gluon distribution function and gluon Sivers 
function. We discuss our choice of WW function and Sivers function in the following subsection. 

C. Sivers function and Weizsacker- Williams function 

In our analysis, we have chosen a kinematical configuration in which proton with momentum P is moving along z 
axis and is transversely polarized in y direction so that 

S ■ (P x k± g ) = k ±gz = cos(j) k± (25) 
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where, k± g = /c ± (cos0 fci , sin0 fci , 0). 

For kj_ g dependence of the unpolarized pdf's, we use a simple factorized and Gaussian form [l2| 

f a /p(x g ,kx g ) = f 9 / P ^ 9 )-^r-^ klJ(kls) - (26) 

We also need a transverse momentum dependent WW function. In our analysis, we have used two choices for it: 
1) A simple Gaussian form as above : 

/ 7/e (x 7 ,fc ±7 ) = / 7/e (3; 7 ) 2 - e -»'-r/<^). (27) 



2) A dipole form 



1 N 

f 7 / e (x 1 ,k± 1 ) = f 7 / e (xj) — — — - 2 (28) 



where iV is a normalization constant, which gets cancelled in the asymmetry. 
For Sivers function we have used two models in our analysis 
• In most part of the analysis we use model I [11] which is given by, 



A N f g/pt (x g ,k ±g ) = A N f g/p ,(x g ) — ^ h(k ±g ) e- fc V< fc i*> cos(&J (29) 



where the gluon Sivers function, A N f g / p t (x g ) is defined as 

,N . 



A N / S/P tfe) = 2Af g (x g )f g/p (x g ) (30) 

where Af g {x g ) is an x-dependent normalization for gluon to be chosen so that the gluon Sivers function obeys 
the positivity bound 

\^ N fg/pt{x g ,k_L g ) 



2 fg/p( x g, k± g ) 

and 



< 1, Vx 9 , fc ±3 (31) 



h(k ±g ) = V2-e'^e~ k ^/ M ? (32) 
Ml 

where the gluon Sivers function, A f g / p t(%g) is given as in model I and Mi is parameter obtained by fitting 
the recent experimental data corresponding to pion and kaon production at HERMES and COMPASS. 

The corresponding parameterizations for quark Sivers function Af u (x) and Afd(x) have been fitted from SIDIS 
data and are given by [l5[ , 

Af f (x) = N f x"f(l - xp { f+ f) . (33) 

a f a fb f f 

where af,bf,Nf for u and d quarks are free parameters obtained by fitting the data. However, there is no 
information available on Af g (x). In our analysis we have used two choices |15| 

(a) N g {x) = {N u {x)+N d {x))/2, 

(b) Af g (x) =Af d (x). 

We also compared the results with model II [3, Gil gi ven by 

A N f g/p ,(x g ,k ±g ) = A N f g/p ,(x g )-^e~ k ^/^ fc f ^ 2 cos(0 fe J, (34) 



and Mq = J(kl g ) where the gluon Sivers function is given as in model I. 

Model I has been used in analysis of SSA in the SIDIS [l2| and DY process es flilj and model II has been used for 
the quark Sivers function to estimate SSA in D meson production at RHIC [la ]. A comparison of k± dependence 
of Sivers function in the two models is given by Fig. 7 of Ref. [IJ] . We will give numerical estimates of asymmetry 
in photoproduction of charmonium using both of these models and will also compare both the parameterizations for 
Af g (x g ) in our analysis in the next section. 
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III. NUMERICAL ESTIMATES FOR THE ASYMMETRY IN J/tp PRODUCTION 

We will now estimate the magnitude of asymmetry using models I and II for both the parameterizations (a) and 
(b) and Gaussian form for WW function k± -dependence. 



The values of best fit parameters of Sivers functions we have used are 31 [ 

N u = 0.40, a u = 0.35, b u = 2.6 , 
N d = -0.97, a d = 0.44, b d = 0.90 , 

Ml = 0.19 GeV 2 . (35) 

These parameters are from new HERMES and COMPASS data [S3, HI] fitted at < Of >= 2AGeV 2 . 

The value of (kj_ g ) is chosen to be the same as (fcj_) for quarks obtained in Ref. [14| by analysis of Cahn effect in 
unpolarized SIDIS from data collected in different energy and Q 2 ranges assu ming a constant Gaussian width. The 
value of (kj_ ) has been chosen to be comparable to (k± g ). As pointed out in [39(, (fc^) in fact depends on energy. 
However, in this paper, we keep this value fixed as above as we do not expect a large variation in this parameter 
related to charmonium production in the model concerned (see below). 

A. Asymmetry using model I 

In model I, the exponential nature of the function h(k±) enables one to perform the k± integration analytically and 
we obtain the expressions for numerator and denominator in asymmetry as 



d 3 ^ d 3 o± 1 



4m 7 



dM^A JV / s / P tK)/ 7/e (x 7 )V2e 



and, 



dyd 2 q T dyd 2 q T s J im 2 Jg/P ' ™"7/«v~7; ~ Ml 

(fci) 2 exp[-4/((fcj) + (fcj 7 ))] ,. 7 ™ m2 , m 

d 3 a^ d 3 ^ 2 ' 4m " 



dM fg/p{x g )f 1 / e (x 1 ) 



dyd 2 q T dyd 2 q T s J im 2 

+ (k 2 ))] 

\M 2 ) (37) 



exp[- g |/((fcj 3 ) + (fcj 7 ))] „ 



where 



(k 2 ) Ml (kl g ) 



B. Asymmetry using Model II 



(38) 



With the k± dependent unpolarized pdf's and gluon Sivers function in Eq. ()26[) and (|34|) numerator becomes 

d 3 aT d 3 a± _ 1 

dy d 2 q T dy d 2 q T s 

x exp[-k± g /a] exp[-fci 9 //3] exp[-g|/a] exp[(2g T k ±g cos((/) k± - 4> qT ))/a] 



f 4m n f k -*-g m 


r-2-K 




/ dM 2 / 


dk± g / a 




)Am 2 c JO 


Jo 





'k ± k ±g A N f g/pt (x g ) f l/e {x 1 ) 



and the denominator becomes, 

I dM / dk± g dcf) k± k ±g f g/p {x g ) / 7 / e (x 7 ) 

4m? JO JO 



dy d 2 q T dy d 2 q T s 

x exp[-kj_ g /a] exp[-kj_ g /P] exp[-g|/a] exp[(2 q T k Xg cos{(/) k± -cj) qT ))/a 



1 



, a ^ cc "(M 2 ) (40) 
7T a p 
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where a = ((q T - k ±g ) 2 ) , (3 = (k\ g ) and M = ^ (k\ g ) = yfp. 

We have used kj_ g ____ = 0.5 GeV, observing that the effective intrinsic motion is limited to k x < yj (k'j_) for 
Gaussian distribution"!?!!] . 

We have estimated the asymmetry with both kinds of parameterizations [labeled (a) and (b)]. The estimates are 
obtained using GRV98LO for gluon distribution functions and Weizsaker- Williams function for photon distribution 
pjil ]. The scale of the gluon distribution has been taken to be §(= M 2 ) (35|. It is worth pointing out that the scale 
evolution of the TMD's including the Sivers function has been worked out in [iO, E| and recently it has been noted 
in [42|, |43[ that in SIDIS the evolution indeed affects the Sivers asymmetry. In fact the evolution effectively produces 
a change in the Gaussian width of the TMDs depending on the scale. The possible effect of such scale dependence is 
not included in our this, exploratory, study of the asymmetry in J ftp production and is deferred to a later publication. 
However, since we are using CEM at LO, the hard scale is between Am 2 and Am 2 D independent of the beam energies 
and the TMD evolution effects within this rather narrow range will be similar for all experiments considered. Our 
results are presented in Figs. (l)-(9). The main features of these plots are summarized below. 

In Figs. (l)-(5) we have shown the asymmetry (A*^ ^ qT ^ s ') as a function of rapidity y and qr respectively for 
JLab (y/s = 4.7 GeV), HERMES (y/s = 7.2 GeV), COMPASS (y/s = 17.33 GeV) and eRHIC (y/s = 31.6 GeV and 
y/s = 158.1 GeV) energies. The plots are for models 1(a) and 1(b) for two choices of the gluon Sivers function. We 
obtain sizable asymmetry in the kinematical regions of all the experiments for model I (b) . The asymmetry is smaller 
in model 1(a). The asymmetry decreases with increase of y for JLab. At HERMES it remains almost independent 
of y for negative rapidity and decreases with increasing y for positive rapidity. For COMPASS as well as for eRHIC 
lower energy the asymmetry increases with y, reaches a maximum and then decreases. This maximum is reached at 
y R3 0.6 for COMPASS and at y ~ 1.2 for eRHIC for model 1(b). For the higher energy at eRHIC no maximum is 
seen in the rapidity region plotted. The asymmetry increases with the increase of y. Both model 1(a) and model 
1(b) show similar qualitative behavior. The asymmetry increases with qx for both models, and for higher values of 
qr ~ 0.6 — 0.7 GeV, it becomes relatively steady. 

In Figs. 6 and 7 we have compared the results for the asymmetry for models I and II at COMPASS energy 
(y/s = 17.33 GeV). Figure 6 shows the asymmetry for models 1(a) and 11(a) and Fig. 7 shows the same for models 
1(b) and 11(b). In model II the asymmetry is larger in magnitude. It increases steadily with qx in model II unlike 
in model I where it first increases rapidly then becomes relatively stable with increase of qr- The difference in the 
qr dependence in model I and model II arises due to the different k± dependence of the Sivers function in the two 
models. Figure 8 shows a comparison of the y and qr dependence of the asymmetry at JLab, HERMES, COMPASS 
and eRHIC. We have used model 1(a) for this comparison. 

So far we calculated the asymmetry using a Gaussian ansatz for the transverse momentum dependence of the WW 
function given by Eq. (|27|). We also do a comparative study taking a different ansatz, namely, the dipole form with 
k = 0.7 GeV. 

We have shown the asymmetry as a function of rapidity or COMPASS energy in Fig. 9 for both Ansatze and using 
models 1(a) and 1(b). It is seen that the asymmetry does not depend much on the choice of the WW function. 

IV. SUMMARY AND CONCLUSION 

We have estimated the magnitude of SSA in electroproduction of J/ip using Weizsacker- Williams equivalent photon 
approximation and models of Sivers function proposed earlier to estimate the asymmetry in SIDIS and Drell-Yan pro- 
cesses. We used the color evaporation model for charmonium production for the numerical estimate of the asymmetry. 
Corresponding to each of the two models of the Sivers function, results have been given for two parameterizations 
of the gluon Sivers function. We have found sizable asymmetry in the energies of COMPASS, HERMES, JLab 
and eRHIC. We have used two different ansatz for the transverse momentum dependence of the WW function and 
found that the asymmetry does not depend much on the choice. Our results based on CEM indicate that it may be 
worthwhile to look at SSA's in charmonium production both from the point of view of comparing different models of 
charmonium production as well as comparing the different models of gluon Sivers function used for estimating SSA 
in other processes. 
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FIG. 1: (Color online) The single spin asymmetry A^ 1 ?T s for the e + — ► e + J/ip + X at JLab as a function of y (left 
panel) and qr (right panel). The plots are for model I with two parameterizations (a) [solid red line] and (b) [dashed blue line]. 
The integration ranges are (0 < qr < 1) GeV and (0 < y < 0.25). The results are given at *Js — 4.7 GeV. 
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FIG. 2: (Color online) The sing le spin asymmetry A°^ 9T for the e + p T -> e + J/ip + X at HERMES as a function of y 
(left panel) and qr (right panel). The plots are for model I with two parameterizations (a) [solid red line] and (b) [dashed blue 
line]. The integration ranges are (0 < qr < 1) GeV and (0 < y < 0.6). The results are given at y/s — 7.2 GeV. 
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FIG. 3: (Color online) The single spin asymmetry A^ Ly,yqT for the e + p r ->• e + J/tp + X at COMPASS as a function of y 
(left panel) and qr (right panel). The plots are for model I with two parameterizations (a) [solid red line] and (b) [dashed blue 
line]. The integration ranges are (0 < qr < 1) GeV and (0 < y < 1). The results are given at yfs = 17.33 GeV. 




FIG. 4: (Color online) The single spin asymmetry A s ^ 1 ^ q ' r ^ for the e + — > e + J/ip + X at eRHIC as a function of y 
(left panel) and qr (right panel). The plots are for model I with two parameterizations (a) [solid red line] and (b) [dashed blue 
line]. The integration ranges are (0 < qr < 1) GeV and (0 < y < 1). The results are given at = 31.6 GeV. 
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FIG. 5: (Color online) The single spin asymmetry y4^ r n ''*' !r for the e + — > e + J/tp + X at eRHIC as a function of y 
(left panel) and qr (right panel). The plots are for model I with two parameterizations (a) [solid red line] and (b) [dashed blue 
line]. The integration ranges are (0 < qr < 1) GeV and (0 < y < 1). The results are given at -Js — 158.1 GeV. 




FIG. 6: (Color online) The single spin asymmetry A*^^' 1 ' 1 ' for the e + — > e + J/ip + X at COMPASS as a function of y 
(left panel) and qr (right panel). The plots are for two models I [solid red line] and II [dashed blue line] with parameterization 
(a). The integration ranges are (0 < qr < 1) GeV and (0 < y < 1). The results are given at y/s = 17.33 GeV. 
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FIG. 7: (Color online) The single spin asymmetry A^ 1 



-4>s) 



for the e + p 1 * -> e + J/ip + X at COMPASS as a function of y 
(left panel) and qr (right panel). The plots are for two models I [solid red line] and II [dashed blue line] with parameterization 
(b). The integration ranges are (0 < qr < 1) GeV and (0 < y < 1). The results are given at yfs = 17.33 GeV. 
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FIG. 8: (Color online) The single spin asymmetry A 



sin(0g 



-<Ps) 



for the e + — > e + J/ip + X as a function of y (left panel) 



and qr (right panel). The plots are for model I with parameterization (a) compared for JLab (yfs = 4.7 GeV) [solid red line], 
HERMES (s/s = 7.2 GeV) [dashed green line], COMPASS {^s = 17.33 GeV) [dotted blue line], eRHIC-1 (^s = 31.6 GeV) 
[long dashed pink line] and eRHIC-2 (t/s = 158.1 GeV) [dot-dashed black line]. 
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FIG. 9: 



(Color online) The single spin asymmetry A*^ ^ tT 



for the e + p' 



1.5 



J/ip + X at COMPASS as a function of y. 



The plots are for model I with two parameterizations (a) [solid red line and dashed green line] and (b) [dotted blue line and long 
dashed pink line] compared for Gaussian and dipole k± dependence of WW function. The results are given at \fs = 17.33 GeV. 



